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Abstract. It is known that Backlund transformations can be used to generate 
stationary axisymmetric solutions of Einstein's vacuum field equations with any 
number of constants. We will use this class of exact solutions to describe the 
exterior vacuum region of numerically calculated neutron stars. Therefore we 
study how an Ernst potential given on the rotation axis and containing an 
arbitrary number of constants can be used to determine the metric everywhere. 
Then we review two methods to determine those constants from a numerically 
calculated solution. Finally, we compare the metric and physical properties of our 
analytic solution with the numerical data and find excellent agreement even for a 
small number of parameters. 



PACS numbers: 04.20-q, 04.25.D-, 04.40.Dg, 95.30.Sf, 97.60.Jd 



1. Introduction 



Solving Einstein's field equations in the vacuum with the assumption of axisymmetry 
and stationarity amounts to little more than solving Ernst's equation, which is a non- 
linear par tial dif f erentia l eq uation. Ernst's formulation of th e stationary axisymmetric 
problem (Ernst (jl968al jbl). Kramer and Neugebauer ( 1968h ) allowed a reformulation 
of the s t ation ary axisymmetric Einstein equations in terms of a "linear syste m" (s ee 
MaisonI (|l978f ). [Belinski and Zakharov (19781 .Harrison UJ28,), NeugebaugB (|l979f) ). 
such that the integrability condition of the latter equals Ernst's equation. 

During the past two decades, several papers have been published focussing on 
constructing analytical Ernst p otentials and applying t hem to various settings (e.g. 
Manko and SibgatuUinl (Il993l) . Stiite and Camenzindl ( 2002 ). Berti and Stergioulad 
( 2004 ) and recently IPappas ( 2009t )). These include solutions to generalizations of 
Ernst's equation, which permit an electromagnetic field and use given axis data in 
the form of rational functions with an ar bitrary number of constants on the axis of 
symmetry based on Sibgatullin's method (jSibg atuUinl (|l99ll) l 

Here, we will present an alternative approach using Backlund tr ansformations 
acting w.l.o.g. on th e Minkowski seed solution (see Neugebauer ( 19791 ). in particular 
Neugebauer ( 1980 ot )) to find concise formulae for the determination of the involved 
parameters. This approach will result in an extremely simple, but in our opinion at 
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the same time elegant recipe, that does not require any sophisticated mathematical 
background. 

This method will then be used to calculate an analytical solution of the field 
equations in vacuum, which approximates the exterior metric for a given model of a 
rotating neutron star very accurately. 

Two different approaches to compute the axis potential are considered. The 
resulting analytical solutions arc then compared with rigorous numerical results using 
the AKM programme (|Ansorg et all (|2003[ l) in the whole vacuum region. We find 
excellent agreement and convergence even for a small number of parameters. 

This paper is organized as follows. In the first part, we review the properties 
of stationary and axisymmetric spacetimes and the relation between the metric and 
Ernst's potential. After that, we introduce the linear problem and show how to solve 
it using Backlund transformations with an arbitrary number of constants, referred 
to as Backlund parameters. Then we provide concise formulae determining those 
parameters from the axis potential. In the next part we illustrate two methods to 
compute the axis potential, namely by using Geroch-Hansen multipole moments and 
alternatively by fitting a quotient of polynomials to given axis values. After that we 
compare our analytical results with full numerical ones. 



2. Construction of the metric by Backlund transformations acting on a 
seed solution 



2.1. The Ernst equation 

The exterior vacuum gravitational field of a stationary axisymmetric uniformly 
rotating matter configuration can be written in Weyl-Lewis-Papapetrou coordinates 

ds' = e-2^ [e^'^ {dg' + dC') + g'dip'] - e^^ {dt + a dipf , (1) 

where U denotes the "Newtonian" gravitati onal potential, a the "gravitom agnetic" 
potential and k the superpotential, cf. e.g. iNeugebauer and Meinell (|2003h . These 
three potentials are functions of g and C alone. The metric ([T]) admits an Abelian group 
of motions with the Killing vectors ^* = SI, ^'^^ —1 asymptotically as + cx) 
and 77' = (5^, 77*772 > for p > 0. 6j is the Kronecker delta, so that has only 
a t-component, and 7/* points in the azimuthal (^-direction. Note that U and a can 
be related to scalar products of the Killing vectors by the coordinate independent 
expressions e^^ = and a = —e~'^rii£^^. 

Ill the vacuum, Einstein's equations for U and a are equal to the so called Ernst 
equation: 

K (/) (^f,,, + fx< + ha) = /,e + fx- (2) 
The metric coefficients a, fc, U and Ernst's potential are related via 

f{gX)^o'^ + ib, 
where b was introduced as 

a^g = ge~^^b^Q, a,(^ = —ge'^^^b^g (3) 
and k can be found by 

k. = Au'g-U', + \e~'''{bl~b',)], (4) 
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It follows from that a,g^ = a^(^g 
obtained by line integration from /. 



and fc k cn are fulfilled, so a and k can be 



2.2. The linear Problem 

It is well-known (see e.g. Neugebauer ( 1980ih . Neugebauer and Kramer ( 1983h ) that 
Ernst's equation is the integrability condition ^^zz — ^,zz of the linear problem 
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Here, we introduced the complex coordinates z :— g -\-\C,^ z := g — iC, and the spectral 
parameter A := {K — if) / {K + iz), where K is an arbitrary complex constant. 
$(A, z, z) is a 2 X 2 matrix, depending on the coordinates as well as on the spectral 
parameter, whereas A and B depend only on the coordinates g and C. One may 
always multiply ([5]) from the right by a so-called gauge matrix which depends only 
on K. Equating all coefficients of A in the integrability condition leads to a set 
of 4 differential equations. One may introduce a potential / in such a way that 
B = f^z/ {f + f) and A — f^^/ (/ + /)• The remaining two equations are then the Ernst 
equation 



2.3. Solving the linear system 

Solving a linear system, which depends on an additional spectral parameter, cor- 
responds to solving a Riemann-Hilbert problem. The idea is then to discuss <i> as a 
holomorphic function of A for fixed z and z and to calculate A and B afterwards. 
With the special ansatz that will be discussed in the following, no integrals have to be 
evaluated. Instead, solving the Riemann-Hilbert problem can be reduced to solving a 
linear system of equations. This special ansatz is to formulate a set of properties, such 
that if a matrix $ has these properties, it automatically fulfils the linear problem. 

A matrix <I> th at has the properties ^■s solution of the linear system 

I 'Neuaebauer i '199d) ). 

= Q + Ai?, = S' -I- T/A, (6) 

where Q, R, S and T are matrices that do not depend on A. 

^ (A) V (-A) 
X(A) -x(-A) 



*(A) 



(7) 



V'(A-i)=x(A) (8) 
^(A = -l,z,z) = x(A = -l,z,z) = l (9) 
For a detailed discussion see e.g. iMauched (|2008[) . iNeugebauer and Krameii dull) 



used the ansatz 



$ = P$o = fiP^o, with /i = 



K 



K 
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to seek new solutions $ from old ones $oj where 

2n 



k=0 



is a finite matrix polynomial in A. We restrict ourselves w.l.o.g. to the Minkowski 
seed solution 



$n = 



Xo -Xo / V 1 
The fundamental lemma of algebra gives 

2n 

detP(A) = det^PiA* = /3( 

i=0 

where the zeros A^ are functions of z and z: 

Afe - 



1 

-1 



fc=i 



Kk - iz 



Kk + iz 
with complex parameters K^. 

bmce det$(Afc) = ^2(ietF(Afc)det<I>o = 0, there is a non-trivial eigenvector for 
each k. Solving this system of algebraic equations using Cramer's rule, one finds the 
n-soliton solution for the particular Minkowski seed solution $o: 
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Here, the Backlund parameters 
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ak = 



X(-Afc)+x(Afc) 

X(-Afe) - x(Afc) 



(11) 



were introduced. Equation (j8]) restricts the complex parameters to Ki = Ki and 
aiO-i — 1 ov Ki = Kk and aiOik — 1- The Ernst potential can be found by setting 
A = 1 in (E 



/(^',C) = x(A = l,^?,C)• (12) 

The solution (|12l) of the Ernst equation was derived by Neugebauer ( 19791) by an 
iterated application of Backlund transformations to the Minkowski space ( Neugebaueil 
Explicit formulae for the metric functions can be found in Appendix 
[Appendix A[ 
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2.4- Recovering the Ernst potential from its axis data 

Thus far, we discussed how to generate general solutions to Ernst's equation in the 
whole (p, C)-plane. It is well-known, that the axis values of the Ernst's potential 



uniqu ely determine the Ernst potential in the whole plane (of. iHauser and Ernst 



The expression for / found in the last section is a quotient of two polynomials 



in Ai, and on the axis, this expression becomes a rational function in In this 
subsection, we will explore the correspondence between the so far (apart from the 
reality condition ([5])) arbitrary Backlund parameters Ki, oci and a given Ernst potential 
on the upper ^-axis ^ denoted as 

f (0 = ^"+^^^-^^^^" . m (13) 

In the following, we will refer to Ernst potentials, that are rational functions on the 
axis as "rational potentials". The formulae that will be present ed in this subse ction 



were found by Neugebauer. They have also been discussed in iMauchen (|2008f ) and 
found recently application in Neugebauer and Hennid ( 2009h . 



From ([5]) together with the definition of A and B as functions of the Ernst 
potential, one finds that the matrix $ factorizes into a matrix depending solely on 
C times a matrix D[K) depending only on K. The matrix D{K) can be evaluated 
explicitly at the branch point W := {\~Q,K — (1,q — 0), where A becomes single 
valued. From PH)) at A = 1 follows that the parameters Ki are the zeros of the 
polynomiai \ [Z (C) ^ (C) + ^ (C) N (C)] . Comb ining these two results, one finds the 
mentioned decomposition into K and C dependent parts on the upper axis A'^: 

$(A^C)= ( ^'j^j \ ^D{K) and 

r,(r^._}_(N{K) + N{K) N{K)~N{k) 
^ ^ " 2i^" Z{K)-Z (K) Z {!<) + Z (K) 

This allows us to show the following: 

Let f (<^) be a rational potential of the form iliA} with arbitrary parameters Oi, hi. 
Let furthermore Ki ^ Kj and Ui ^ Uj Vi ^ j- For any set of parameters {oi, bi\ there 
is exactly one set of parameters {Ki^Ui], such that the Backlund transformation B 
acting on the Minkowski seed solution results in an Ernst potential f 

B[fo^l,K,,a,]^ f 

that becomes the given rational potential on the axis. The Backlund parameters can 
be explicitly calculated by evaluating the expressions 

- [Z (C) N (C) + Z (C) N (C)] - n - ' (14) 

1=1 



Proof: 
Evaluation of 



N {K,) Z (K,) 



<i>(x,)(i>o)-M f 



I One has to consider the complex continuation of the real functions Z(Q and N(Q. 
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leads to (|15p . This overdetermined system of equations for the has solutions, iff 
Z (Ki) N (Ki) + Z [Ki) N [Ki] = 0. Another way of proving this fact is by simply 
evaluating the expressions (|15p and (jl4p using the representation of / in form of 
determinants (e.g. (ITUl) ) at A = 1, and verifying that they are fulfilled. 
q.e.d. 

The general algorithm to determine / (p, C) from f {g = 0, C) works as follows: 
First of all, one has to calculate the Ki by evaluating The fundamental lemma 

of algebra guarantees existence and uniqueness. Then one can easily find the by 
evaluating (jTS]) . The equations and (jTS]) are universal in the sense that one does 
not need to know beforehand whether the Ki or the ai are complex or real, nor does 
the order of the polynomial matter. Equation (I15[) allows the calculation of even 
if Ki is a zero of N or of Z . The result can be easily generalized by using Bernoulli's 
rule for coinciding zeros {Ki = Kj for i ^ j). The have to be treated as usual 
functions of Ki. 



3. Different approaches of constructing rational Ernst potentials on the 
rotation axis 



3.1. Reflectional symmetric Ernst potentials 



In Newtonian theory one can s how that rota t ing fl uid bodies in equilibrium always 
have reflectional symmetry (cf. iLichtenstein (Il933l) ). The same symmetry is to be 
expected in Einsteinian theory (cf. iLindblomI ( 19921 ) i and we assume it to exist for our 
purpose. Therefore we study the i nfluence of reflectional symmetr y bef ore constructin g 
rational potentials. As shown in iMeinel and Neugebauer (Il995[) and iKor dai (Il995l) . 
this kind of symmetry of the metric is uniquely characterized by a simple relation for 
the Ernst potential / on the upper part A'^ of the C-axis: 

(16) 



/(C)/(-C) = 1- 



Condition (ITCl) is fulfi l led fo r aU C, iff N{C) = Z{~C). As shown in the proof of 
Theorem 2 of iKordasI this is equivalent to the following conditions for the 

coefficients: 



, _ , _ , I ai, n even 

On = —an, On-l = On-l . . . Ol — < _ , , 

I — ai, n odd 

So we have reduced the 2n complex parameters to 2n real ones. 



(17) 



3.2. Construction from multipole moments 



As a first method of generat i ng a r ational Ernst potential /(C) on the rotation axis, 
we refer to iManko and Ruiz ( 1998h . They start from the expansion of the function 
X{Q at C = oo: 



1 - /(C) 
1 + /(C) 



(18) 



where the coefficients nij determine the Geroch -Hansen mul t ipole moments Pj (the 
explicit relations between rUj and Pj are given in lFodor et al.1 ( 19891 )). After plugging 
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in the ansatz (1131) and evaluating the algebraic set of equations, which results from 
comparing the coefficients of the powers of up to n = 2N in (jlSp . one ends up with 
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(19) 



and an analogue formula for N{() by changing the sign "— " to "+" in the first column 
of the numerator. From this formula one can see that condition (|17p is automatically 
fulfilled when TOj is real for even j and purely ima ginary otherwis e, which is another 
well-known formulation of reflection symmetry (see iKordas ( 1995( )). 

Since only far-field information is employed by use of the multipole moments, it 
is guaranteed that the /(C) constructed in this way will coincide well with the solution 
from which the multipole moments were taken in a neighbourhood of infinity. 



3.3. Fitting a given (non-rational) Ernst potential 

Another method is constructing a rational potential /(C) from a given function g((^) 
on the rotation axis, which is especially useful in a numerical context. There, points 
Q on the axis are selected where /(C) shall exactly coincide with g{C). If the points 
Ci are suitably chosen, we find a good approximation / of g on the whole axis. The 
same ansatz (1131) is employed. 

A major advantage of this method is that /(C) and g{() will also coincide near 
the star's surface, if one or some of the Q's are situated there. 

Rational potentials found this way do not automatically have the expected 
behaviour at large distances (all multipole moments and especially mass M = too = Po 
and angular momentum J — —imi — —iPi of the analytical solution may vary from 
the given numerical ones). To guarantee the right asymptotic behaviour, we treated 
the coefficients in (jl3|) as functions of m„ at fixed M and J. The to„ were varied until 
/(C) and g{() coincided at the gridpoints Ci- 



4. Results 



After constructing a rational Ernst potential one has to calculate the Backlund 
parameters according to formulae (|14|) and (|15|) . For polynomials in numerator and 
denominator of degree n there are 2n real parameters. The analytical approximate 
solution (often referred to as Backlund solution, e.g. as presented in formula (|10l) ) is 
a 2n-fold Backlund transformation acting on the Minkowski seed solution. 

In the following we present solutions with n = 1, 2, 3 and 4 to show the increasing 
accuracy of the approximation when taking more and more parameters into account. 



For our numerical calculations we used the AKM programme (jAnsorg et al 



where different coordinates had been introduced. We have chosen to compare 
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our results in Weyl-Lewis-Papaetrou coordinates and have taken the necessary 
transformations into account. 



4.I. Comparison of metric functions 



As a first resuh we have calculated the metric functions and the relative difference 
between analytical and numerical values (see figure [1]). For the numerical calculations 
one has to choose an equation of state. In this section we decided to model our neutron 
star to be composed of matter of constant energy density. The star we have chosen 
to treat as an example has a ratio of polar to equatorial radius of 0.7 (this ratio is 
always given in the coordinates used in Ansorg et al.l (|2003f ) ) and a norma lized central 
energy density of 1 in order to compare the results with Nozawa et al.l (l998il and 
Ansorg et all (|2003l ). 



• n = 1 



■AKM 



0.80 



1.00 




1.20 



0.55 



0.30 



0.05 



<l 

O 




3 

(a) Function e'^^ and relative 
difference for g = • ■ ■ 5g, 



3 

(b) Function e^* and relative 
difference for g = Qc ■ ■ ■ 5ge 



3 

(c) Function a and relative 
difference for q = Qc ■ ■ ■ 5ge 



Figure 1: Metric potentials as functions of coordinate radius q in the equatorial plane 
(homogeneous equation of state and fitted rational potential on the axis with fixed M and 
J). 



These approximate solutions were constructed by fitting the given numerical Ernst 
potential on the axis, holding M and J fixed. For n = 1 (i.e. a twofold Backlund 
transformation) this already determines the two available real parameters so that no 
additional gridpoint can be choosen. This is equivalent to approximating the neutron 
star's exterior metric by a Kerr solution with the same mass and angular momentum. 
We expect therefore the approximation to coincide well with the numerically calculated 
potentials only near infinity, as we already explained in the end of section 13. 2| and to 
differ most near the star's surface. 

In a 2n-fold Backlund transformation with n > 1 one has to fix more parameters. 
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In fitting the Ernst potential on the axis, our algorithm chooses a new gridpoint where 
- in the previously calculated (2n — l)-fold transformation - the approximation and the 
numerical potential differ most. For a fourfold Backlund transformation, this is the 
star's surface, while for higher transformations the additional gridpoints are usually 
well spaced on the axis, which shows the efficiency of our algorithm. 

As one may observe e.g. in figure [1] one gains about one order of magnitude in 
accuracy by going from a n-fold to a (n + l)-fold Backlund transformation, i.e. with 
two parameters more. 



^.2. Comparison of physical properties 

As we approximate the exterior vacuum region of a numerically calculated spacetime 
it is necessary to define what physical properties mean in our context. Of course one 
can simply read off mass M and angular momentum J from the far-field behaviour of 
the metric, as well as higher multipoles. This is sufficient if one only wants to analyse 
far-field behaviour. 



■71 = 1 



n = 2 



n = 3 



■• n = 4 



■AKM 



0.138 



0.131 

CD 

0.125 



0.118 



CD 

<l-4 
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^-6 



7r/4 



7r/2 



7r/4 



7r/2 



Figure 2: Potential e^^ and relative difference as a function of i9 := arctan(£i/("). e^^ 
coincides at pole and equator (for n > 1), because we obtained the analytical vacuum metric 
by fitting the numerical Ernst potential at the rotation axis (and the first gridpoint is the 
pole i9 = 0). 



A more intriguing question is how to determine the star's surface in the appro- 
ximated solution. Using the co-rotating potential V = U' defined by (see e.g. 
Meinel et al.1 (|200^ ) 

2^^' - [(1 + na)^ - n^e^e-^''] , 



e 



the star's boundary is given hy V — Vq — const. However, given only the exterior 
metric (and not the rotation speed fl of the star), it is not possible to construct the co- 
rotating potentials directly. Our approach to deal with this problem is the following: 
Since on the rotation axis g = the co-rotating potentials coincide with the non- 
rotating ones, we read off the approximation of Vq from the approximated potential U 
at the north pole Vq = U'{g = 0, C = Cp) — U{q = 0, ^ = Cp) a-^d use the constancy 
of Vo to calculate the rotation speed 57 at the equator 

e^^o =e2^'(^. = £<e,C = 0). (20) 



The values Qc and Q-p are supplied by the numerical solution. Once is known we 
can evaluate V = U' in the whole vacuum region and V" = Vo is valid at the pole and 
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equator by construction. Of course, the so constructed V will not be constant over 
the original (numerical) star's surface, and also Vq will differ from the value in the 
numerically calculated solution if U does not coincide at the pole. This can be seen 
in the change of behaviour from n = 1 to n = 2, since for n > 1 a fitting gridpoint is 
located at the star's pole. 

To illustrate this behaviour we plotted the potential e^^ along the star's surface 
and the relative difference between the 2n-fold Backlund solution and the AKM 
solution in figure [21 

To remedy also this non-constantness, we simply define the star's surface in our 
approximation as the region where e^^ = e^^° . One can see in figure [3] that this area 
is almost identical with the star's surface for solutions generated by 2n-fold Backlund 
transformations with n > 1. This way of interpreting the analytical solution has the 
advantage that one only needs two surface points (the pole Cp and the equator g^) and 
does not need to worry about the whole surface function gs{C)- 



71 = 1 n = 2 n==3 n = 4 AKM 



1.00 




0.00 I 1 1 ' 

0.00 0.33 0.67 1.00 

Figure 3: Area e^^' = e^^ [(1 + Qa f - Q,'^ g'^e"'^^] = e^^o in comparison with the star's 
surface. 



Finally we calculated the fo llo wing physical quanti ties in table [T] (formulae can 
be found in lNozawa et al.l (|l998[ ) or lAnsorg et all |2003h ): 



-Rcirc equatorial circumferential radius 
Zr, polar redshift 

equatorial redshift in the forward direction 
equatorial redshift in the backward direction. 

Analogous results for different equations of state can be found in Appendix 
[Appendix B[ 

4-3. Critical points and singularities 

A closer look at the spectral parameter A = y/ (K — iz) / (K + iz), where the 
analogously defined are part of the formulae for the metric functions, indicates 
some critical points, because the square root has a branch cut which we took along 
the negative real axis. One easily sees that the points 

Zfci := iKk = + i^Kk, 

where the radicand of Afe vanishes and 
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Table 1: Detailed comparison of physical quantities (homogeneous equation of state 
(/J, = fio) and fitted rational potential on the axis with fixed M and J). Here, pc ~ Pc/moj ^ ~ 

= ^?circ A*o ^'i'^ ~ "^A'o are normalized values of the physical 
quantities. The column labelled AKM shows numerically calculated values, whereas BT[rz] 
refers to values of a 2n-fold Backlund solution. ABT[n] labels the relative difference between 
BT[n] and AKM. 



Model 


AKM 


BT[1] 


BT[2] 


BT[3] 


BT[4] 


ABT[1] 


ABT[2] 


ABT[3] 


ABT[4] 


Pc 


f.OO 


















rp/rc 


0.700 


















n 


1.4f 


1.58 


1.42 


1.41 


1.41 


1.2e-01 


7.4e-03 


1.4e-03 


2.9e-04 


M 


0.136 


0.136 


0.136 


0.136 


0.136 










J 


0.0141 


0.0141 


0.0141 


0.0141 


0.0141 












0.345 


0.337 


0.344 


0.345 


0.345 


2.4e-02 


4.6e-03 


9.3e-04 


1.9e-04 




1.71 


1.91 


1.71 


1.71 


1.71 


1.2e-01 








^cq 


-0.163 


-0.295 


-0.170 


-0.164 


-0.163 


8.1e-01 


4.3e-02 


8.5e-03 


1.8e-03 


^oq 


11.4 


11.7 


11.1 


11.3 


11.3 


2.6e-02 


2.8e-02 


5.7e-03 


l.le-03 



where the radicand of Afe goes to infinity can become problematic. For example 
in the case of a 4- fold Backlund transformation with 4 complex Ki {Ki, K2 — 
Ki, if 3 = —Ki and Ka^ = —Ki) the metric function e^^ has a jump along 

As the formulae for the calculation of the metric functions contain determinants 
in the denominator, even singularities of the analytic solution may occur. Therefore 
we have numerically calculated the zeros of the denominator (see figure|4]) and roughly 
speaking they usually appear near or between the Ki. 




0.00 0.33 0.67 1.00 0.00 0.33 0.67 1.00 

Q/ Se el Qc 

Figure 4: Star's surface (sohd hne), critical points (triangles) and singularities (boxes) of 
the Backlund solutions. The exterior solution is well behaved, because the difficult points 
are inside the star and the singularities are situated near the critical points. 

All these problems are controllable for our purpose. As we search for an analytical 
vacuum solution for the exterior of a neutron star, this is well behaved as long as the 
critical points and singularities are located in the interior region of the star. To decide 
whether our solution is applicable or not, we simply added a plot of the star's surface 
with the critical points to the output of our approximation programme. 
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4.4- Numerical methods 



To ca lculate the approximation, we augmented the AKM programme (lAnsorg et al 



to incorporate the formulae given above. Two main tasks had to be 



done: Calculation of determinants (e.g. (jA.l[) . (jA.2p and (|A.3p ) and root finding for 
polynomials p4|) . 

To calculate the determinants, Gaussian forward elimination was used to convert 
the matrix to triangular form, whence the determinant is given by multiplication of 
the diagonal elements. For our purposes, this method was found to be numerically 
stable by using partial pivoting, i.e. transposing columns until the largest element in 
each row was on the diagonal. 

Polynomial root finding was accomplished for n = 2 and n = 3 by direct solving 
the quadratic or cubic equation, while for n > 3 we determined the roots as the 
eigenvalues of the companion matrix. However, Gaussian elimination was not sufficient 
in this case and we had to resort to a QR deco mposition of the companion matrix 



(adapted from the GSL (cf. iGalassi et ail (|2009l ))) 



At last, to find the singularities, we had to search for zeros of the denominator 
in (jA.ll) . Since the determinant viewed as function in the complex (p, C)-plane, 
has branch cuts where it is not even continous, the obvious Newton's metho d fails 
spectacularly. The Nelder-Mead simplex algorithm (jNelder and MeadI (Il965l) ') does 



work for non-continuous functions and is well behaved for this problem. We start 
with a number of random trial points in the {p, ^)-plane, from which this algorithm 
searches local minima of the absolute value of the determinant. Afterwards we check 
if these minima really are zeros. 



4-5. Conclusion 

It is remarkable that the whole exterior region of a stationary and axisymmetric 
neutron star can be expressed with an accuracy greater than 99% everywhere using 
only 6 to 10 real parameters (pc and (^p to determine the star's surface and 4 to 8 
multipole moments, m„ or values of Ernst's potential on the rotation axis). As there 
are many ways to construct a rational potential on the axis one may choose whether 
near or far-field information should be considered. From this point of view, there is 
no "best approximation", but rather it can be adjusted for specific needs. 

In each case one has to be aware of the singularities and as long as they are 
located in the interior of the star, there is no problem. 



Appendix A. Explicit formulae for the determination of metric functions 



An alternative formula (obtainable from pop by performing the limit A 1, some 
manipulation of the determinants and inserting the definition of r^) for the calculation 
of Ernst's potential is 



/(g,C) = --, AT = 



1 







K" 



Oiiri 
^a2r2 



K. 



2n 



K. 



n-l 



,(A.l) 
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with r,; = \i{Ki + \z). Of course, is the real part of / and one could obtain a 
from its imaginary part b via ([3]). A straightforward calculation leads to the explicit 
formula 



(a-ao)e2^ 2n£» + 25ft I — 



—nz —Q—[n—l)z — i 
iff airiK'^-^ Kl'-^ 



(A.2) 



K. 



2n 



a2r2K^' 



K" 



where the integration is already done. The only thing left to be done is the 
determination of the constant ao, which can be calculated by demanding that a — > 
as g ^ oo. 

Furthermore the last metric function e^*"' could be calculated via line integration of 
([4]) . By reformulating th e linear problem ([5]) and some manipulation of determinants 
(cf. E ramer et alJ l|l991 ). [Krameil (|l980 )) one ends up with the following formula: 

(A.3) 



^2{k-k„) _ 



n 



rj{Kj +iz) 



PIM2, 



(Ki + iz)" 
{K2 + iz)" 



{K2n + '^zr a2nr2nK^~ 



a2r2 



and D2 is the same as Di with replaced by l/a^. Finally the integration constant 
ko can be evaluated by demanding that e^*^ — 1 as gi — > cxd. 



Appendix B. Results for different equations of state and approximation 
methods 

In order to avoid a flood of data we do not present the results for three different 
equations of state (strange quark, constant energy density and polytropic) with three 
different approximation methods (via moments m„, fitting a nonrational potential and 
fitting a nonrational potential with fixed M and J). The following figures and tables 
provide one example of each method, but not necessarily the best one for the given 
equation of state. 



Appendix B.l. Polytropic stars approximated by moments m„ 

As an example with polytropic equation of state (polytropic exponent F = 2, 
polytropic constant K) we have chosen p,c = 1 and r^/rp, = 0.834 in ord er to compare 
the va lues for the physical quantities with Nozawa et al.l (fl998.1 and Ansorg et al. 
1 20031) . 

A first look at figure IBll shows, that all approximated functions coincide with 
errors less than 1%. The eightfold Backlund solution is missing from those data. 
Since critical points and singularities appeared outside the star, we had to discard this 
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n=l • ■•n = 2 n = 3 n = 4 AKM 




12345 12345 12345 

g/ Qc el Qc qI Qc 

(a) Function e^^ and relative (b) Function e^* and relative (c) Function a and relative 
difference for q = Qc . . . 5gc- difference for g = Qc ■ ■ ■ 5gc- difference for g = ■ • • 5pe- 

Figure Bl: Metric potentials as functions of coordinate radius q in the equatorial plane 
(polytropic equation of state and rational potential via moments m„ on the axis). 



solution. As these approximated solutions are generated only by far-field information 
we observed quite generally that critical points appear more often outside the star. 

As this approximation method does not include any information about the star's 
surface, at first sight one could think that e^^ is not approximated very well. The 
graphs (see figure IB2[) however show quite a good approximation even in this case 
(especially for n > 1). Remembering the definition of Vq and fl (cf. (HH])) it is clear 
that the values of e^^ for -i? = and d ~ tt/2 coincide. 



n = 2 



n — 3 



• n = 4 



• AKM 
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Figure B2: Potential e^^ as a function of ■& :— arctan(f(/("). 



Figure IB3I confirms that the real star's surface coincides well with the one 
determined in the approximation by e^ 
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n=l n:^2 n = 4 AKM 
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Figure B3: Area e^^' = e^^ [(1 + Qaf - Q.'^g^e''^"] = e'^° in comparison with the star's 
surface. 

Finally, we calculated some physical quantities in table IBll 

Table Bl: Detailed comparison of physical quantities (polytropic equation of state 
= p + ^p/K) and rational potential via m„). Here, jlc = fJ^cK, Q = QK^^^, M = 
M/K^^^ , -Rcirc = Rciic/ K^^^ and J = J/K are normahzed values of the physical quantities. 
The column labelled AKM shows numerically calculated values, whereas BT[n] refers to 
values of a 2n-fold Backlund solution. ABT[n] labels the relative difference between BT[7i] 
and AKM. 



Model 


AKM 


BT[1] 


BT[2] 


BT[3] 


ABT[1] 


ABT[2] 


ABT[3] 


Ac 


1.00 














rp/re 


0.834 














n 


0.400 


0.425 


0.400 


0.401 


6.1e-02 


1.8e-03 


2.4e-04 


M 


0.161 


0.161 


0.161 


0.161 








J 


9.49e-03 


9.49e-03 


9.49e-03 


9.49e-03 








-Rcirc 


0.679 


0.678 


0.679 


0.679 


2.1e-03 


8.2e-05 


5.8e-06 




0.458 


0.467 


0.458 


0.458 


2.0e-02 


l.Oe-03 


8.0e-05 




-0.0601 


-0.0845 


-0.0594 


-0.0602 


4.1e-01 


1.2e-02 


1.6e-03 




1.04 


1.09 


1.04 


1.04 


4.4e-02 


1.8e-03 


1.8e-04 



Appendix B.2. Strange quark stars approximated by fitting Ernst's potential at the 
rotation axis 

As an example with strange quark eqtaation of state (linear relation between energy 
density and pressure, MIT-Bag constant B) we have chosen Pc = 2 and rp/r-p = .5 in 



order to compare the values for the physical quantities with lAnsorg et aL ( 20031 ). As 



the star is quite relativistic, figure [B4l shows that all approximated functions do not 
reach the accuracy of the other examples. 

The next figure shows the area e^^ = e^^". The fact that the curve for the n = 1 
Backlund solution does not end at the equatorial radius, might be confusing at the 
first sight (especially when in figure [B5] the error of e^^ vanishes there). As the surface 
area e^^ is not necessarily a monotonic function, there can be two intersection curves 
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Figure B4: Metric potentials as functions of coordinate radius q in the equatorial plane 
(strange quark equation of state and fitted rational potential on the axis). 
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Figure B5: Potential e^^ as a function of ' 



arctan(g/(^). e^^ coincides at pole and 



equator, because we obtained the analytical vacuum metric by fitting the numerical Ernst's 
potential at the rotation axis (and the first gridpoint is the pole i9 = 0). 



that need not be connected. In figure [B6l this is the case, where one intersection curve 
ends at g ^ 0.67 Qc and the other starts at g — Qq. 

Finally, we calculated some physical quantities in table IB2I 
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n = 1 n = 2 n = 3 n = 4 AKM 
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Figure B6: Area e^^' = e^^ [(1 + f2o)^ - n^p^e"*''] = e^^" in comparison with the star's 
surface. 

Table B2: Detailed comparison of physical quantities (strange matter equation of state 
(/i = 3p+4B) and fitted rational potential on the axis). Here, pc = Pc/B, Q = M = 

MB^^"^, RciTc = i?circ B^^^ and J = JB are normalized values of the physical quantities. The 

column labelled AKM shows numerically calculated values, whereas BT[n] refers to values of 
a 2n-fold Bficklund solution. ABT[n] labels the relative diffcrcucc between BT[jil and AKM. 




Mo(l(>l AKM BT|i| BT|2| BT|:i| BT|1| ABT|il ABT|2| ABT|:i| ABT|ll 



Pc 2.00 
rp/re 0.500 





3.43 


4.20 


3.57 


3.48 


3.45 


2.2e-01 


4.0e-02 


1.5e-02 


4.3e-03 


M 


0.0355 


0.0293 


0.0349 


0.0354 


0.0355 


1.8e-01 


1.6e-02 


3.6C-03 


8.3e-05 


J 


l.lOe-03 


7.12e-04 


1.06e-03 


1.92e-03 


l.lOe-03 


3.5e-01 


4.1e-02 


9.6e-03 


3.3e-04 


-Rcirc 


0.141 


0.129 


0.139 


0.140 


0.141 


8.4e-02 


1.9e-02 


7.8e-03 


2.4e-03 




0.726 


0.726 


0.726 


0.726 


0.726 










7f 


-0.302 


-0.411 
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-0.312 
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3. 6c- 01 
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